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1. 


Introduction. 


The  finite  element  method  has  become  the  main  tool  in  computational 
mechanics.  The  MAKABASE ■ [ 1 ] ,  [2]  contains  about  20,000  references  on  finite 
element  and  2000  boundary  element  technology.  Recently  the  new  direction  in 
the  finite  element  theory  and  practice  appeared,  the  p  and  h-p  versions, 
which  utilize  high  degree  of  elements.  About  3-4  dozen  references  about 
p  and  h-p  versions  are  available,  all  of  them  related  to  the  elliptic  prob¬ 
lems.  For  the  survey  of  today’s  state  of  the  art,  we  refer  to  e.g.  [3]. 

The  present  paper  addresses  the  basic  problems  of  the  p-version  for  the 
parabolic  equation  with  both  variables,  x  and  t  discreted  via  p-version. 
It  concentrates  on  the  case  when  in  the  time  variables  only  one  interval  is 
used.  The  paper  gives  the  error  estimates  and  presents  some  numerical 

*V  .'t.' 

aspects.  We  restrict  themselves  to  the  basic  features  of  the  method.  Various 
generalizations  will  be  presented  in  forthcoming  papers. 


2.  The  p-version  for  the  initial  value  problem  for  an  ordinary  differential 
equation 

2.1.  Preliminaries  and  problem  formulation. 

Let  I  =  (-1,1),  I  =  [-1,1],  t  e  I,  X  »  Lgd)  be  the  usual  space 
furnished  with  the  norm 

1 

(2.1)  ||u||x  = 

Let 

£’  =  {v  e  C°°U)  i  v(l)  =  0), 

where  c'dl)  is  the  usual  space  of  functions  with  all  continuous  derivatives 


u2dt 


l/Z 


on  I.  For  any  A  >  0  and  v  e  we  define 


(2.2) 


x*Xv"x  • 


•  dv  A) 

where  we  denoted  v  =  Let  be  the  completion  of  C  with  respect  to 


the  norm  II  •  | 


V 


Lemma  2.1.  Let  v  e  and 


(2.3) 


IMI*  =  +  l|Av,lx- 


Then 

(2.4) 


<  ||v||  <V2||v|| 

A  A  A 


Proof.  We  have 


imi; 


’-i 


)  dt  =  [(*]  +  ^v^2]dt  “2 


^  +  Av |  dt  = 


,1 


-1 


vvdt 

-1 


=  imi: 


Y(v2(t))dt  = 

at 


-1 


=  llvll?  +  v2(-l)  >  ||V||2  . 
A  A 


On  the  other  hand 


imi: 


s  2("s«x*»x'''x)  * 


2iivii: 


Lemma  2.1  implies  that  the  spaces  Y  and  Z..  are  equivalent. 

A  A 


On  XxY  ,  u  e  X,  v  e  Y  we  define  the  bilinear  form 


V  W  W  TV  •)!  V ■  Jl  V-  ■J.'TV  ■.  'VT^^T^*  J^VT  <jv  ■/■  '.-V  i-vjv  <jv  <;vm 


(2.5) 


B  (u.v)  = 


f  u(- 1  *  Xu)dt' 

J-1 


Further,  let  F  e  Y'  be  a  linear  functional  on  Y  . 

A  A 


We  will  define  now 


Problem  T  .  Given  F  e  Y^,  find  6  X  such  that 


B^CUq.v)  =  F(v)  V  v  e  Y^. 


Theorem  2.2.  Problem  has  a  unique  solution  un  €  X  and  ||Ug||x  <  ||F||V/. 


Proof.  By  Schwarz  inequality 


(2.6a) 


|Ba(u,v)|  <  ||u||xl|v||Y  . 

X 


Given  u  e  X,  then  there  is  v  e  Y^  such  that 


A  +  XV  =  «. 


Obviously  ||v||  =  ||u||  .  This  yields  immediately 

YA  x 


(2.6b) 


inf  sup  |B. (u, v) |  >  1 
ueX  veYx  A 


I U|| X=1  II v|| Y  <1 


and  analogously 
(2.6c) 


inf  sup  | B  (u, v) |  >  1 . 
ueX  A 


veY. 


I  v  ||  Y  =1  II  ull  X<1 


Now  Theorem  2.2  follows  immediately  from  Theorem  5.2.1  of  [4], 


Let  us  show  that  the  solution  of  the  problem  T  is  a  weak  solution  u 

A  0 


of  the  initial  value  problem 


to® 


> 


w 

% 

‘\A 

Si 

V, 

v, 


,yj 

» 


?<■ 


>3 

4 


ft 


s 

y-j 


■$? 

Si 


$ 


l 


3 


(2.7a) 


(2.7b) 


A"  +  Au0  =  f 


uQ  ( - 1 )  =  aA 


F(v)  =  fvdt  +  av( -1 ) . 


(2.7c) 


To  prove  it,  it  suffices  to  show  that  if  f  €  C^(I)  and  e  C^(I! 

(2.7a,b)  then  uQ  is  the  solution  of  the  problem  9  . 


solves 


We  indeed  have  for  any  v  e 


Wv>  - 


fr*S) 


vdt  +Uq(-1)v(-1)-  =  fv  dt  +  av(-l). 


This  shows  that  any  classical  solution  (i.e.,  a  solution  belonging  to 

C1(T) )  is  also  the  solution  of  the  problem  9^  with  F  given  by  (2.7c)  and 

hence  the  problem  9^  is  a  weak  formulation  of  (2.7a,b). 

We  mention  that  because  of  Lemma  2. 1  and  the  Sobolev  imbedding  theorem 
0  -  r* 

v  €  c — »C  (I)  and  hence  F(v)  =  J  fvdt+av(-l)  e  Y^.  F(v)  can  be 

obviously  identified  with  a  function  from  H  *(I).  We  note  that  u  e  X  has 
in  general  no  trace  in  t  =  -1,  but  we  still  have  an  initial  value  problem. 
The  initial  condition  is  now  of  the  type  of  the  natural  boundary  condition  and 
not  the  essential  one. 

2.2.  The  p-version  of  the  finite  element  method  for  the  problem  9  . 

A 


(2.8) 


Let  q  >  1,  integer 
)  S  =  Sq_1  =  {u 


€  X  |  u  is  a  polynomial  of  degree  q-  1} 


V  =  §)q  =  (v  €  Ya  |  v  is  a  polynomial  of  degree  q) . 


_  V  —  •*  f-  ^  V  <■  ^  *r  rn  rn  *"(.  if".  m  *"*»  ■,  •'V  »’V  ' 


Theorem  2.3. 


i)  Let  u  e  S,  v  e  V,  then 


(2. 10) 


I  ( u ,  v )  [  —  11  *~M x 11  '  ”  Y 


ii) 

(2. 11) 


1  - 1/2 

d^(q)  =  inf  sup  |B^(u,v)|  > 

lies'1-1 

l|v|lY  =1  luljjSl 


iii)  Let  u  e  S,  u  *  0,  then 


(2. 12) 


sup  | B  (u, v) |  >  0. 
=  1 


Proof. 


i)  (2.10)  follows  from  (2.6a). 

ii)  Let  P  ,  k  =  0,1,...  be  the  Legendre  polynomials.  Then  v  e  V  can 
be  written  in  the  form 

q 

(2.13a)  »  *  Z Vk 

k=0 


with  the  constraint 

(2. 13b) 


q 


k=0 


Let  the  X-orthogonal  projection  operator  of  X  onto  S  and  let 

2  =  Viw’ 


w 


V 

—  +  A  v 
A  ’ 


v  €  V. 


Then  for 


UmiBillVIfJ  VLW>  ?M  I'J.'^yT^JTVVvyv  W  'A  -  V.’.V^.V  V,  V  V  V 


VVVTWWU’TOTXTT'TT 


Hence 


q 


q-1  20? 


q_1  20? 


=  y  run  <  fi »_ai_i  y  i<2(y^ 

/L  2k+l  ~  [  2 (2q+l ) J  zL  2k+l  "  Zl  2k+l 


k=0 


k=0 


k=0 


and  therefore 


Ba,z,v)  -  I|£||2*A2 


2 £2 

,. v  2  .2  V  k  „ v  2  ^  1  ....  ,2 

Z.vdt  =  ll-IL  +  A  >  ?rr— r  ^  II -II „  +  —  II Av|| 


'1 


AX 


2k+l  ~  AX  2q 


-1 


k=0 


“  lq^V^Y  ^by  Lemma  2-  1  ■ 

A 


Realizing  that  ||z||  =  B  (z,  v),  (2.11)  follows. 

X  A  - 

iii)  Given  u  €  Sq  *,  let  v(t)  =  f  udt.  Then  v  e  §)q  and 

Jt 


(2. 14) 


and  (2.12)  follows. 


B  (u,v)  = 


,1 


-1 


r  2 


Avvjdt  >  i||u||2  >  0, 


We  can  now  define  the  p-version  of  the  finite  element  method  for  the 
problem  T  :  Given  F  e  Y'  and  q  >  1,  integer,  find  u  e  1  such  that 

A  A  CJ 


(2. 15) 


B,  ( u  ,  v )  =  F(v)  V  v  €  §)q  =  V. 
A  q 


Theorem  2.3,  together  with  Theorem  6.2.1  of  [4]  yields 

Theorem  2.4.  There  is  a  unique  u^  satisfying  (2.15).  If  u^  e  X  is  the 

exact  solution  of  the  problem  f  ,  then 

A 


(2. 16) 


lu0“uJx  -  H+2q1/2]  inf  || uQ  -  w|| x_ 
q  wgS 


The  p-version  of  the  finite  element  method  reduces  to  the  solution  of  a 
system  of  linear  equations.  Let  P^(t)  be  the  Legendre  polynomials,  then  we 


can  write 


q-1 

u  » y\ 

q  zL  k 


v  ■  Xvwv *  v- 


Hence  (2.15)  reduces  to  a  system  of  linear  equations  for  a^: 
q-1 

Z“AlPk'pj.rpj)  -  F(pj.rpj>’  J  - 0 . o- 


Of  course,  we  can  use  another  basis  function  of  S  and  V,  too. 

-I/O 

If  (2.11)  is  optimal,  i.e.,  d(q)  <  Cq  ,  then  (see  [5])  there 

exists  u  e  X,  which  is  solution  of  the  problem  9  for  certain  F  e  Y' 

A  A 

such  that  the  finite  element  solution  diverges,  i.e. 


u-u  „ — as  q — >  co . 
q  X  ^ 


On  the  other  ha^d,  for  a  restrictive  class  y  of  the  solutions  u  of  the 


problem  9  ,  we  have 
A 


|U-  u  I! x  <  C  inf  II u  -  w|| x 
q  weS 


with  C  independent  of  q  or  growing  with  q  much  more  slowly  than  q 
In  the  next  sections  we  will  address  more  of  these  aspects. 

2.3.  The  optimality  of  (2.11). 


q-1 

u  * 


WWVJ re 


f>  n  I 


© 


£  £  p/kc 

j=0  k=0  J  K  J-1  J 


pq-i  ?«2 

^.2  V  ^  2 

Zl  2k+l  2q+l 


(2. 17) 


q  q-1 


B  (u.v)  = 


£  E  Vi  f  q-l  ? 

j=0  k=0  *  J  J-1  J  k  V  2  jf 

;  +  xz. 


Obviously 


R  =  A  sup 


^  2a  .0. 

,y  .  j  j| 

'Z-  2J+1 

j=o 


J  (2J+1) 


I  2a“/(2j+l)=l 


(2. 18) 


inf  A 


q-l 

2  V  02  2 

L  P0  (2j+l) 

J=o 


where  inf  is  taken  over  8.  such  that 

J 


(2. 19) 


Select 


f^'1  282  f^'1 

>2  y  y 

ZL  20  +  1  2q+ 1  Z. 

\J=0  \}=0 


*0  =  l. 


(2. 20) 


I3J  =  ^(Sj+l)  0  =  0 . q-l 


.  i  A.  — N»-  s-.”V  A»  l  L.Pm.  t  /*j«  i.* »  u* j|  _*j 


*r  »,  %  .V.v^v.’v  yv 


with  K  determined  so  that  (2.  19)  holds.  Then 


9-2  9 

C^q  <  K  ^  <  C2qJ, 


with  C^,C2  independent  of  q  and 

Q  S  Cq_1, 

where  C  is  independent  of  q.  Assuming  that  (2. 19)  holds,  it  is  obvious 
that  there  is  a  function  0  <  ft^(x)  <  oo,  0  <  x  <  oo  such  that 

llvlly  SJMq^+l. 

A  A 

Further  there  is  a  function  0  <  Jf2(x)  <  w,  0<x<oo  such  that 


sup  |B^(u, v) |  <  Jf2 ( q)  —  +  R. 

“j 

2 

£  2a,/(2j+l)«l 


Hence 

.  ,  ,  .  _  ID  ,  ,  ,1  _l/2  - 1/2 . ,  ^l(q),-l/2 

d  (q)  =  inf  sup  |B  (u,v)|  <  (tf  (q)-  +  C  q  )(1+ - =-) 

II  V||  Y  =1  Hull  X<1  *  A2 

and  we  have  proved 

Theorem  2.5.  If  A  >  C_1/2q1/2K2(q) ,  A  >  J<j/2(q),  then 

.  ,  ,  „  ..,1/2  -1/2 

d^(q)  <  4C  q 

where  C  is  the  same  constant  as  in  (2.20),  independent  of  q. 

We  see  that  for  A  which  is  large  in  comparison  with  q  the  estimate 
(2.11)  is  optimal,  i.e., 


1 im  sup  d  (q)  <  Cq 
A->oo  A 


It  does  not  mean  that  for  fixed  A  and  q,  d  (q)  is  small.  In  fact  we  can 

A 

expect  that  d  (q)  achieves  minimum  for  some  finite  q  depending  on  A. 

A  U 

The  constant  d  ( q )  can  be  computed  numerically.  Let  u  €  S,  v  e  V, 

A 


(2.21) 


u<t)  =  Z“jVi<tK 


(2.22) 


=  T.0 


where  p1  =  p^U)  =  J  j  =  2,3 . q,  and  P^U),  j  =  0,1... 

are  normalized  Legendre  polynomials  of  degree  j.  Then 


where 


(2.23) 


B^(u, v)  =  0  Aa, 


ai  bl 

C2  a2  b2  ° 


0  -  ’  q-1 

c  a 

q  q 


i  =  1,2 . q 


[2i-l ) (2i+l 


i  =  1,2 . q-1 


v.v.v.v.v.v;.h'.-.v.v,-,vv,.* 


v(2i-l ) (2i+3) 


“l . “q] 


*1 . 8,1 


T 

=  a  a 


T 

=  3  Cg, 


el  fl 


®2  f2 
d3  e3 


fq-4  eq-3  dq-2  Sq-2  fq-2 
fq-3  ®q-2  Vl  ®q-l 
fq-2  eq-l  dq 


3+6A2+4A2 


1  ,  2A 

X  (2i  +  l) (2i-3)  ’  1 


3,  i=2 . q-1 


'(2i+3)(2i-l)(2i+l) 


.  .  .  ,q-2. 


V.V.V.V.V'.V.AV.V/AY.V 


/./.v.v.v.v .v.v 


■  »“-V- -«r  rv 


rn'tVIii'l'iIVl'i 


.•iMk’Jk'A1  .1 


T  T 

,  (s  _  .  „  IB  AA  |3 

dA(q)  lnf  .  T  T  ,1/2 
/3TC|3<1  (/3  AA 


-  .nr  |*W2 

(3  L 


and  hence  d^(q)  is  the  smallest  eigenvalue  of  the  problem 

AAT/3  =  d^(q)C0. 

Table  2.  1  presents  the  values  of  d^(q)  for  various  q  and  A.  Figures 

2.1  and  2.2  present  d  ( q )  in  the  dependence  on  q  and  A.  We  see  that  for 

any  fixed  A,  d  (q)  first  decreases  with  q  and  then  increases,  for  fixed 
A 


q,  d^(q)  decreases  with  A  and 


..  V'1.1 

Jil  W 


The  detailed  theoretical  analysis  of  the  structure  of  d  (q)  is  not 

A 

avai lable. 
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. 98680 

.79701 

.71727 

.67610 

.66514 

.66381 

. 66237 

. 66203 

.66159 

5 

. 99506 

. 82899 

. 67937 

. 58675 

. 56138 

.55829 

. 55493 

.55415 

.55312 

10 

. 99875 

. 92492 

. 74648 

. 51874 

. 44260 

. 43325 

. 42311 

. 42076 

. 41764 

15 

. 99944 

.96118 

. 82747 

. 53586 

.39710 

.37941 

. 36024 

. 35582 

. 34994 

20 

. 99969 

. 97678 

. 8784F 

. 58089 

.38175 

.35411 

. 32406 

.31715 

. 30796 

25 

. 99980 

. 98468 

.91095 

. 62875 

.38291 

. 34423 

.30172 

. 29194 

. 27897 

30 

. 99986 

.98918 

.93257 

. 67097 

. 39403 

.34403 

.28781 

. 27482 

. 25763 

40 

. 99992 

.99381 

.95813 

. 73776 

.43085 

.36091 

. 27550 

.25519 

. 22834 

50 

. 99995 

.99601 

.97179 

. 78735 

. 47194 

. 38901 

. 27593 

. 24743 

.20947 

100 

. 99999 

. 99899 

. 99240 

.91037 

.62570 

. 52393 

. 34373 

.28026 

.  17611 

Table  2.1.  The  values  of  d  (q). 

A 


2.4.  The  set  of  perfect  solutions. 

We  have  shown  in  section  2.2  and  2.3  that  there  are  solutions  u  e  X 


such  that 


(2.25) 


llu-u J 

.  ,  „  q  *  =  R  .(u) 

inf  ||u— w||x  q,X 

weS 


can  be  arbitrarily  large  provided  that  q  and  X  are  sufficiently  large.  On 
the  other  hand,  we  have  proven  that 

1/2 

R  .(u)  <  Cq1'* 
q,  A 


where  C  is  independent  of  q  and  u.  Let  ^p(x)  <  »,  1  <  x  <  oo  be  a  non- 
—  ■ - K  -.0  as  x — > oo.  Then  if  D  >  0  is  arbi- 


decreasing  function  such  that 


trary,  there  exists  u  e  X  such  that 

sup  R  (u)#>  *  (q)  >  D. 
q,A  q’ 

Hence  we  can  ask  to  characterize  the  set  r(<p)  c  X  depending  on  <p  such  that 

sup  sup  R  ,(u)f  *(q)  <  oo. 
uey  q, A  q’ 

The  set  Tf(<p)  will  be  called  the  set  of  ^-perfect  solutions.  An  especially 
important  case  is  y(^)  for  <p  =  1.  It  is  not  easy  to  give  a  precise  charac¬ 
terization  of  the  set  y((p)  although  its  importance  is  obvious.  We  will  give 
only  some  sufficient  conditions  and  based  upon  numerical  experiments  we  will 


also  formulate  a  conjecture. 


u  €  X,  u  =  ^<x 


.P..  Then 
J  J 


where 


inf  || u  —  w ||  x  =  || u  -  u || x» 
weS 


*  *>  ’ 


Let 


u 

q 


be  the  ignite  element  solution  defined  by  (2.15)  and 


u  -  u 

q 


u  -  u  -  z 


with 


z  =  u  -  u  e  S. 

q 


Then 


B  (z,v)  =  B.  (u  -u+u-u,v)  =  B.(u-u.v)  V  v  e  V. 
A  A  q  A 

q 

For  v  . 

k=0 


B^(z,v)  = 

Because  ( by  Lemma  2.1) 

A2^  <  A2||v||2  <  ||v||2 

\ 

we  get 

IBx(z,v)|  <  |aq|||v||Y 

and  hence  using  (2.11) 


(u-u)^-  ^+Avj 


dt  =  A/3  a  . 

q  q 


1/2 


II z ||  <  2  | a  |q 


Iu-Vx s  (fi*j 

v  1=rr 


oo  .1/2  r  .  .  1/2 

— r  ?  I  Ct  I  q 

>  oc  +  2 1 a  | qiX^  =  inf  ||u-w||  1+2  —3 - — 

—  J  <*  J  u*=5  X  f  "  si 1/2 


which  leads  to  the  following 


Theorem  2.6.  Let  <p(x)  be  a  nondecreasing  function 


R) 

lJ=q 


,  U  €  X,  U  =  y^Ot. 


,  P,  and 
k  k 


(2.26) 


r  |  1  1/2  . 

f  l<Xqlq  -1 

sup  1+2  — f - Y75  V  (q)  <  “• 

<?  .  r;  si 1/2 


(i  «5) 

LJ=q  JJ 


u  €  y(<p). 


Consider  some  concrete  examples: 


“i  -v^-ip0m-4  £ 


pt(t) 

/§k5T(2k-l)(2k+3) 


and  hence  u^  is  a  ^-perfect  solution  for  ip  =  1. 


ii)  Let  0  <  p  <  1  then 


_L_  =  f  i  p*. 


In  this  case  Theorem  2.6  gives  no  indication  whether  u^  is  a  <p-perfect 

1/2 

solution  for  any  <p  (with  exception  of  course  of  <p  =  x  ). 

The  observation  about  perfect  solutions  is  practically  important. 

Theorem  2.6  shows  that  if  the  solution  is  unsmooth,  the  finite  element 

solution  has  essentially  the  same  accuracy  as  the  best  approximation  also  for 

large  q  which  is  needed  to  get  an  acceptable  accuracy.  If  the  solution  is 

smooth  then  the  acceptable  accuracy  is  achieved  for  small  q  and  hence  the 
1/2 

factor  q  is  not  important. 


v.v/f.V 


X."  O-  *  "  < 


2.5.  A  numerical  example. 


Let  us  consider  the  problem 


(2.27a) 


A  +  Xu  =  * 


(2.27b) 


u(-l)  =  0. 


The  solution  u  is 


u  ( t )  =  1  -  e 


-X2(t+1) 


Assume  that  the  basis  functions  of  S  and  V  are  given  by  (2.21)  and 


(2.22),  respectively.  Then  the  p-version  reduces  to  the  solution  of  linear 


system 


Act  =  F 


where  A  is  given  by  (2.23)  and 


F  =  [A,-  —.0.  .  .0)T. 


The  values  of  the  function  R  (u)  given  by  (2.25)  are  given  in  the  Table 

q,  A 


2.2.  Computations  were  performed  in  double  precision.  Figure  2.3  presents 


R  ..(u)  in  dependence  on  the  polynomial  degree  q  for  various  A.  The  shown 

q.  a 


slope  is  the  theoretical  one  based  on  Theorem  2.4  (n  =  0.5). 


1 


i 

0 


1.01317 


1.00493 


10  1.00000 


15  1.00000 


20  1.00000 


25  1.00000 


30  1.00000 


40  1.00000 


50  1.00000 


100  1.00000 


3 

5 

10 

20 

25 

1. 14719 

1. 11079 

1 . 03955 

1.01089 

1.00705 

1. 14853 

1. 19126 

1.09510 

1.02909 

1.01910 

1.07101 

1.20361 

1.23124 

1. 10008 

1.06915 

1.03765 

1. 14908 

1.28929 

1. 18268 

1. 13488 

1.00018 

1. 10947 

1.28816 

1. 25360 

1.20130 

1 . 00000 

1.08241 

1.26821 

1. 30273 

1. 25826 

1 . 00000 

1 . 02000 

1.24639 

1. 33032 

1. 30113 

1 . 00000 

1 . 00000 

1.20763 

1. 34308 

1. 34560 

1 . 00000 

1 . 00000 

1. 17530 

1. 33333 

1.35521 

1 . 00000 

1 . 00000 

1 . 00000 

1.26922 

1. 31420 

1.00279  1.00179  1.00045 


1.00767  1.00494  1.00125 


1.02947  1.01926  1.00495 


1.06241  1.04165  1.01103 


1. 10258  1. 07027  1 1. 01934 


1. 14601  1. 1030811.02970 


1. 18923  1. 13803  1.04190 


1.2653711.2073711.07083 


1.32000  1.26778  1.10408 


1.37144  1.3806511.27263 


Table  2.2.  The  values  of  R  ,(u)  for  various  q  and 


LINE  OF  MAX’S 


Figure  2.3. 


2  3  4  5  6  7890  20  30  40  50 

DEGREE  q 

R  ..(u)  in  dependence  on  the  polynomial  degree  q,  slope 

<4 »  A 

M  =  0.5  based  on  (2.16). 


iwWKiHSrSSWi 


Table  2.2  and  Figure  2.3  show  that  for  large  A,  R  .  (u)  grows,  but 

q,  A 

1/2 

stays  bounded  (realize  that  we  have  used  q  =  100  leading  to  q  =10  and 

dioo^lOO)  =  0.17611  in  Table  2.1).  It  is  interesting  to  observe  that  for 

fixed  A,  R  ,(u)  first  increases  and  then  decreases  to  1.  Hence  Table  2.2 
q,  A 

and  Figure  2.3  suggest 

2 

Conjecture  A.  Let  u^  =  e  *  (t  +  1).  Then  u^  e  y(.<p)  with  <f>  =  1  or 
<p(x)  =  log(x).  In  fact  we  have  seen  that  for  any  A  and  any  q  in  a 
practical  range 

R  v  ( u )  <  1.5. 
q,  A 


2.6.  The  error  estimates. 

We  have  shown  in  previous  sections  that 


II  u 


Uq II X  *  C(q) 


inf  ||  u  -  w|| 
weS 


where  X  =  L2(I).  Assuming  that 

00 

u  ’  Z\Pk 

k=0 

we  have 

1/2 

2 
k 


00 

|u-wlx=  Y, 


and  hence  we  deal  with  the  error  in  L-2 ( I )  of  the  reminder  of  extension  in 
Legendre  polynomials.  There  are  many  results  which  can  be  used  here.  For 
more  see  also  [6].  For  example,  using  the  standard  results  about  Legendre  and 
Jacobi  polynomials  we  have 


M  p  2  ' 

inf  || u -  wH  <  C(p)q"P  (l-t2)pf— 1  dt 

weS  X  J  ,  WpJ 

“  X  J 


where  an  explicit  form  for  C(p)  can  be  given. 

„  _  .  -x2(t  +  n 

u  =  1  -  e  we  get 


For  example  for 


(2.28)  inf  1| u  w|| Y  <  C.CpJq'PxP'1. 

weS  X  1 

If  u  is  an  analytic  function  on  I  then  we  get  exponential  rate  of  conver¬ 
gence.  Let  us  derive  another  estimate  we  mentioned  above.  Using  Taylor’s 
formula  we  get 


(2.29) 


II  u  -  w||x 


2q 


1_ 
q!  ‘ 


We  see  that  for  large  q  or  small  X  the  rate  is  of  order  (q! ) 


2.7.  Additional  aspects. 

We  discussed  the  finite  element  method  for  the  case  when  S  c  X,  i.e., 
we  were  interested  in  the  error  measured  in  the  L^norm.  Because  we  can 
exchange  the  spaces  S  and  V  and  still  keeping  the  same  properties  of  the 
bilinear  form  present,  we  can  seek  the  finite  element  solution  of  the  problem 
(2.7)  also  in  the  form 

q 

u(t)  =  V  a.  P,  (t) 
k=0 

u  ( —  1 )  =  aA 


and  determine  the  coefficients  a^. 

In  our  analysis  we  have  assumed  that  the  solution  u  e  X  and  it  has  no 
trace  and  hence  we  cannot  say  anything  about  u(+l).  We  can  modify  our 
approach  so  that  value  of  u(l)  will  be  also  the  solution.  To  this  end  let 

X[U  =  XxR1,  u  =  (u.p),  u  e  X,  p  e  !R^ 


with  the  norm 


and  let 


be  the  space  without  the  condition  at  t  =  1  with  the  norm 


i 


and  the  bilinear  form  will  be 


B(u, p; v) 


pv(l) 

A  I 

( 

I 


which  leads  to  a  weak  solution  of  the  problem  2.7  with  p  =  u(l). 

We  discussed  only  the  case  of  single  interval.  It  is  possible  to  treat 
also  in  an  analogous  way  the  case  when  the  time  interval  is  derived  into  sub- 
intervals,  etc.  We  will  elaborate  on  these  or  similar  problems  in  the  future. 


3.  The  p-version  for  the  parabolic  problem 

3.1.  Preliminaries  and  problem  formulation. 

2 

Let  n  c  R  be  a  bounded,  Lipschitz  domain  with  a  piecewise  analytic 
boundary  T.  Let  further  D  =  IxQ,  I  =  (-1,1).  Then  we  will  consider  the 
problem 


(3.  1) 


-  Au  =  f  on  D 
o  t. 


u  =  0  on  Ixf 


(3.2)  u(-l,x)  =  g(x)  on  n. 

Let  X  =  X(D)  =  L2 ( I , H 1 ( n ) ) ,  where  H^(Q)  is  the  standard  Sobolev 

space,  as  well  as  H°(Q)  =  L2(fi)  with 


(3.3) 


lull,  =  |  Vu  |  dx 

H!(n)  Jn 


iiuir,  dt. 
A1(n) 


By  H  i£2)  we  denote  the  usual  Sobolev  space  with  the  norm 


|  uvdxl 

Jn 


lull  =  sup  „ - . 

H  (n»  vfi'tO)  #*(0) 


By  6*  we  denote  the  space 

£*  =  {v  e  C°°(  D)  |  v(t,x)  has  for  any  tel  compact  support 
i  n  £2  and  v(  1 ,  x)  =  0} . 

We  denote  by  Y  =  Y(D)  the  completion  of  £’  in  the  norm 


(IMI2,  +  IM|2  )dt. 

h  Mn)  h  (n) 


where  v  =  — .  We  have 
ot 


Lemma  3.1.  Let  v  e  and  V(t,x)  be  such  that  for  any  tel,  V  e  A  (£2) 


(3.5) 


VV(t,x)Vzdx  =  vzdx  V  z  e  A  (£2). 


Ml  =  ||  v|| 

H  (£2)  H1  (£2) 


Proof.  Obviously  V(t,x)  is  uniquely  determined  by  (3.5).  We  have  then 


1  -1 
h  (n) 


I  vzdx| 


o  1  "‘•"ol 

zeH  (fi)  H  (fi) 


I  VVVzdx I 


ol  11  "ol 

zeH  (n)  h  (n) 


H1  (fi) 


Hence  we  have 


(3.6) 


(IIVII^  +l|v||2  )dt, 
H1 ( fi)  fi1 (Q) 


where  V  is  defined  by  (3.5). 

Let  now  (A.,u.),  u  6  H^(fi)  be  the  eigenpair  of  the  eigenvalue  problem 
J  J  J 


-Au  .  =  A  .u  .  on  fi 
J  J  J 


u .  =  0  on  T. 
J 


Vu  .Vzdx  =  A^  u  .zdx  V  z  e  H^fi) 
J  J  J 


and  any  u  e  H  (fi)  can  be  written  in  the  form 


(3.7) 


u  =  >  oc.u. 

L-,  i  i 


r  00 

2  ,  V  2 

u  dx  =  )ai 


|  Vu | ^dx 


Z2.2 

Vi' 


Hence  X  -  X(D)  is  the  set  of  the  functions  of  the  form 


u(t,x)  =  y^a. (t)u. (x] 


(3.8)  II u|| x  =  /  <  +00- 

**-1  i=l 

Further,  it  is  easy  to  see  that  if  v  €  £’  then 


(3.9a) 

CD 

v(t,x)  =  (t)ut  (x) 

i=l 

with 

(3.9b) 

R.(l)  =  0 

and 

IM)2  =  Y/S2(t)^-dt 

Hence 

(3. 10) 

H  vll  Y  =  £(|32(t)^  +  02(t)A2 

J-lli=l  Ai 

and  Y  co  s  of  the  set  of  functions  of  the  form  (3.9a,b)  so  that  (3.10) 
is  finite. 

Remark  3.1.  Because  R^(t)  e  ( I ) ,  i  =  1,2,...,  fL(l)  exists  for 
i  =  1,2 . 

Using  now  Lemma  2.1,  we  see  that  the  norms 

(3. 11) 

and 

(3. 12) 


,>  .1 1.1 


are  equivalent. 


Let  us  define  on  XxY  the  bilinear  form 


(3. 13)  B(u, v)  = 


( -uv+7uVv)dxdt 


1  00 

=  >  (-a.  3.  +  A 2cc.I3. 

Lu  ii  l  i  l 

J-14=l 


.13.  )  dt 


(3.8)  and  (3.10)  yield  immediately  that 


(3. 14) 


|B(u,  v)  |  <  ||u||xl|v||Y 


and  using  the  same  argument  as  in  section  2.1  we  get 


(3.  15a) 


(3. 15b) 


inf  sup  |B(u,v)|  >  d  >  0 

u<=X  ve  Y 

lull  X=1  ||V|IY<1 


inf  sup  |B(u, v)|  >  d  >  0. 
veY  ueX  1 


|v||Y=l  ||u||  X<1 


Let  us  now  define  problem  ?.  Given  F  €  Y'  find  uQ  €  X  such  hat 


(3. 16) 


B(uQ,v)  =  F(v)  V  v  6  Y. 


From  (3.14)  and  (3.15)  we  have  that  problem  f  has  unique  solution  for 


any  F  €  Y'  and 


|uQ|lx  <  CI|F||y, 


where  C  does  not  depend  on  F. 


It  is  also  easy  to  show  that  the  solution  u^  of  the  problem  ?  is  a 
weak  solution  of  (3.1)  and  (3.2)  with 


F(v)  =  g(x)v(-l, x)dx +  f (t,x)v(t,x)dxdt 


in  the  sense  that  if  u^  is  smooth  and  satisfies  (3.1)  and  (3.2)  then  it 
satisfies  also  (3.16). 


K« 


26 


Remark  3.2.  We  restricted  ourselves  to  the  special  case  of  Q  c  IR  with 
Lipschitz  domain.  We  did  it  only  for  simplicity.  For  example  in  this  section 
no  restrictive  assumption  on  £2  has  been  made  at  all. 

3.2.  The  semidiscrete  problem.  Discretization  in  x. 

In  this  section  we  will  consider  the  semidiscrete  (in  t)  solution  of 
the  problem  f.  Let  R  c  ft*(£2)  be  a  finite  dimensional  subspace  of  functions 


S  =  (u  e  X  I  u(t,x)  e  R  V  t  e  I> 


V  =  (v  e  Y  I  v(t,x)  e  R  V  t  e  1} 


Let  u  e  S  be  such  that 
s 


(3. 17) 


B(u  ,v)  =  F(v)  V  v  e  V. 

s 


Then  we  will  call  us  the  semidiscrete  solution  of  the  problem  ?. 


In  (3.4) 


we  have  defined  the  norm  ||  •  | 


Let  u  e  R,  then  we  define 


H  (Q) 


|  uvdx | 


(3. 18) 


Obviously 


(3. 19) 


m  i,2 

11  U|1  -1  =  SUP  IR - ' 

Hr  «!)  vsR  »fil(n) 


Hull  <  Hull 

(n)  h  x(n) 


We  will  say  that  the  space  R  has  property  X  if  there  is  a  number  H(R)  <  1 


such  that 


(3.20) 


|u||  >  H(R)flu|| 

hr  (n)  h  (n) 


holds  for  any  u  e  R. 


Obviously  H(R)  depends  on  R.  We  will  discuss  it  later. 


then  the  norms  IMIY  and  II'IL  are  equivalent. 

R  yr 

Theorem  3.2.  Let  R  have  property  K,  u  €  S,  v  €  V.  Then 


(3.24) 


(3.25) 


(3.26) 


I  B(u,  v)  |  <  ||u||xl|v||Y, 


inf  sup  |  B(  u,  v)  |  >  CpJf(R) , 
ueS  vsV  ■' 

lull  X=1  ||  vj|  Y<1 


sup  |B(u,v)|  >0  V  v  e  V,  v  *  0, 
ueS 

|u|lx=1 


where  and  C 2  are  independent  of  R. 


Proof.  (3.24)  follows  from  (3.14)  and  (3.22).  Let 


VC  V  *  »  / 

u  e  S,  u  =  >  a.  u„  .  , 

Z-.  i  R,  i 


v<  v  *  *  / 

V,  v  =  /  j3.u0  ..  Then 

ZL  i  R.  i 


r1  f^(R)  i  1 

:u'vi  =  X  ViM'  cr*Visi) ' 

-1M  =  1  ’ 


Using  the  results  of  section  2,  given  u^  e  S  there  exists  v  €  V  such 


B(u,v)  =  ||u|i: 


and  using  equivalency  of  the  norms  Y  and  Y  we  get 

K  R 


Hull  "  "Y 

-  NXM7  ranr  a  C"U#**IR) 

R 


s 


which  leads  immediately  to  (3.25).  (3.26)  follows  also  easily  as  in 

Section  2. 

Using  now  Theorem  6.2.1  of  [4]  we  get  immediately 

Theorem  3.3.  Let  R  have  property  K,  and  u^  are  the  solutions  of 
(3.16)  and  (3.17),  respectively.  Then 


|uo-us"x  £  HIrT  ln£  lluo-"'x 

weS 


and  C  is  independent  of  R. 


Let  us  now  discuss  the  accuracy  of  us_  We  have 


B( u  -u  , v)  =  0  V  v  e  V. 
0  s 


Write  u  =  p  +  £,  p,£  €  S,  where  p  will  be  judiciously  selected.  Then  we 


B(£,v)  =  B(uQ-p,v)  V  v  e  V 


B(Uq~p,v)  =  (-(Uq~p)v + V(uQ-p)Vv)dxdt. 


Let  us  select  p(x,t)  =  PgUgCx.t),  where  is  L^-orthogonal  projec- 

O  ^ 

tion  of  H  (fi)  onto  R  (for  every  fixed  t).  Then  for  any  v  e  V 


B(uQ-p,v)  = 


r1  r 

(7(VVo 

-rn 


)Vv)dxdt. 


The  bilinear  form  B(u,v)  satisfies  the  conditions  (3.24),  (3.25),  and  (3.26) 

without  the  property  X  if  using  the  norm  YD  instead  of  Y.  Because 

K 


(7(Vpouo)Vv)dxdt  -  i|u0-P0u0llxl|v|lY 


r  » 


.■.V.VA.V’A'A'.V.V.V.VV.V 


re  ge 


■«x  -  c||uo-pouo"x' 


where  C  is  independent  of  R.  Hence  we  have 

Theorem  3.4.  Let  be  the  solution  of  the  problem  ?.  Then 

■Us-Vx  £  cluo-pouo"x' 


where  by  PQ  we  denoted  the  L^-orthogonal  projection  of  X  onto  R  without 
assuming  the  property  X.  □ 


3.3.  The  condition  X  and  related  results. 

In  the  previous  section  we  mentioned  condition  X  of  the  subspace  R. 
We  will  now  elaborate  more  on  it.  Let 


(3.27) 


Hull 


£(R)  =  sup 


H^n) 


ueR  liu||L2(£l) 


<  +00 


(3.28) 


T)(R)  =  sup  II U  -  Pj^ull L  (n)  <  +co 

uefi^n)  2 

Hull  <1 

ft^D) 


(3.29) 


A(R)  =  sup  II u  P qu|| L  (£2)  <  +oo. 
ueft^Q)  2 

Hull  -  <1 

bhn) 


7),  and  X  depend  on  R.  Here  we  denoted  P^,  resp.  P^,  the  L^-, 
resp.  H1-projection  operator  of  ^(fl)  on  R.  We  will  call  now  R  to  be 
(?.  T),  X)-regular. 


Theorem  3.5.  Assume  that  R  is  (£,  tj,  A)-regular.  Then  R  has  property  X 
and 


H(R)  >  (1+Cn) 


Proof.  Let  u  e  R.  Then 


|  uvdx| 
J  O 


MTj 

H  <n)  vsH'tn)  h1 (n) 


|  uvdx | 

J  n 


HR  V6R  "  V(0) 


Hence ,  we  have 


I  f  uP  vdx| 

Jn  1 

Hull  1  =  sup  — 77— jj - +  sup 

h  (n)  ai,_,  11 

veH  (£2)  H  (n)  veH  ( Q) 


<  ||u||  +  t)(R)||u|| 

HR1(n) 


I  [  u( v— P  v ) dx  | 
Jn  A 


H1  (Q) 


|J  uvdx  |  II  u||  L 


iiu|1  -1  "  SUP  gfRlTvH - 

veR  ?(R)|1  llL0(n) 


Hence 


H  (Q) 

i 

hr  (n) 


S  1  +  $1J. 


Theorem  3.6.  Assume  that  R  is  (£,  7),  A)-regular.  Then  for  any  u  e  H  (Q), 


|u-Pnu[|  <  (1+tj(R)+A(R))?(R). 

0  H^n) 


Proof.  We  can  write 


u  -  P_u  =  u  -  P.  u  +  P.  u  -  P_u. 


Hence 


l|u-P_u||  <  ||u-P.u||  +  (v(R)+*(R))£(R)  *  (1+t?(R)+A(R))£(F0. 

0  ftX(n)  1  ft!(£2) 

Remark  3.3.  In  the  above  theorem  we  do  not  need  to  use  operator  P  .  It  is 
enough  to  use  any  other  one  P^  with  t?(R)  instead  to  tj(R).  Then  Theorem 
3.5  holds  when  t)(R)  is  replaced  by  t)(R).  For  example  we  can  use  some  type 
of  an  interpolation  operator  which  in  the  case  of  unsmooth  domain  Q  yields 
t)(R)  <  T)(R)  . 

The  conditions  (3.27)  -  (3.29)  are  closely  related  to  the  inverse  proper 
ty  and  the  duality  principle  in  the  theory  of  finite  element  method.  Let  us 
mention  here  some  known  results. 

i)  Consider  finite  element  method  with  a  quasi-uniform  mesh.  Then 

?(R)  S  h_1C(p), 

where  h  is  the  size  of  the  element  and  C  depends  on  p  and  the  mesh  but 
is  independent  of  h. 

In  the  case  of  the  p-version  we  have 

?(R)  *  p2Q 

where  Q  depends  on  the  mesh. 

ii)  Assuming  that  £1  is  smooth,  or  all  its  angles  are  <  n,  we  have 

U(R)  *  C(p)h 
or 

tj(R)  <  p  1Q. 

iii)  Finally  without  any  additional  assumption  on  the  mesh  we  get 


A(R)  <  Cp  *. 


Remark  3.4.  In  this  section  we  did  not  use  any  assumptions  about  fi  and  its 
dimension,  except  in  the  last  part,  where  we  listed  the  results  about  the 
finite  element  space  R. 

Let  us  now  analyze  the  p-version  for  the  one-dimensional  (in  x)  case. 
Assume  now  that 

O  J 

R(p)  =  {u  €  H  (I)  |  u  is  a  polynomial  of  degree  p} . 

Then  we  have  (see  e.g.  [6]) 

£(R)  <  Cp2 

2 

and  the  factor  p  is  the  best  possible.  Further  we  have 


||u-P  v||  <  Cp  || v|| 

1  L2  1  H1 ( I ) 

Hence  we  have 


tj(R)  <  Cp" 


(3.30) 


H(R)  >  Cp 


Similarly  we  get 


(3.31) 


'U-P0u|lol 

H 1  ( I ) 


S  Cp||u| 


ftX(  I ) 


Let  us  show  now  that  (3.31)  can  be  improved. 


Lemma  3.7.  Let  s  >  0  integer,  u  €  ft  (I)  and 


f  <u‘ 


s+1 ] ,2.  2 ,s  . 

)  ( 1  — x  )  dx  =  A  <  +oo. 


(3.32) 


[s+1]  dS  *u 


|u-P  u||  <  C(s)p~S+(1/2)A1/2. 

°  fl1(I) 


Proof .  We  can  assume  that  u  is  smooth  because  of  density  of  smooth 


functions.  Assume  that 


-I 


a  (p  x,-P  1  ) 

m  m+1  m-1 


where  Pm  are  as  before  Legendre  polynomials.  Then 


u'  =  )  a  (2m+l)P 
/  .  m  m 


r  * 

P.u  =  Yb  (P  -p  , ). 
0  m  m+1  m-1 


For  the  sake  of  simplicity  we  will  restrict  ourselves  to  the  case  when  u  is 

symmetric.  The  antisymmetric  case  can  be  treated  similarly.  In  our  symmetric 

case  we  have  b  =  a  =  0  for  all  m  even.  Let  m  =  2k-l,  k  =  l,2,...,k_, 
mm  0 

k„  =  £  p  even.  Let  us  assume  that  b  =  a  +  z  .  Because  P„  is 
0  2  m  m  m  0 

L2~projectlon  we  have 

f1  r  k° 

XZ2k-l(P2k_P2k-2}  "  a2k +l(P2k +2_P2k  ]  (Pi+l~Pi-l)dx  =  0 
J-14c=l  u  u  u  J 

i=l , 2, • • • , p-1 . 

This  system  leads  to  the  following  system  of  linear  equations: 

(2+I)zi"|Z3  =  0  (i  =  !) 


cws 


v’lsKvl 


Ki 

K 

B 


v*  ^  O 


2  ,2  2,  2  „  ,. 

"  g +  (5+g)z3  "  gz5  =  0  (l  =  3) 


2  2  2  2 
4k-3  Z2(k+1 )-l  +  ^4k-3+4k+"l^Z2k+l  ”  4k+TZ2(k+l )-l  =  0  (l  =  2k_1) 


2  2  2  2 

4k0-3Z2kQ-3  +  (4k^3+4k^TT}z2k0-l  =  ”  4k0+la2kQ+l  (l  =  2kcf1}' 


This  system  can  be  solved  explicitly  to  obtain 


Z2k-1  =  -Ck0k(2k-1)a2k0+l>  k  =  1  ’  '  '  ’  k0  ’  2k0+1  =  P+1 


in  which 


;  =  4V3 

'k0  8k3.6^-5k0-3’ 


Hence  we  have 


A)  ” 

|u'  *  <P0U,-|L  (I)  *  Z^k-l'4"-”*  I 


k=k0+1 


4-i(4k-') s  40.iko*“' 


where  Q  =  ^  a2k-l^k-^  anc*  ^  *s  bounded  from  above  and  below  indepen- 

k=ko+1 

dently  of  k^.  Because  of  our  assumption  we  have 


(u' )2dx  =  Y  a2(2m+l) 

/  .  m 


Hence 


V a2(2m+l)1+2s  <  C  (u[s+1] )2(l-x2)Sdx  <  CA. 
/  .  m 

m=l  -1 


U2V,)2  S  CA(2k0»l)-1-2s 


>  V  V  V  V  V  V  V  >  *»  *V 


O  v' 


«v 


I  4-, 


(4k-l)  <  Ckr 


k=k0+1 


Therefore 


|u'-,poul'll?2<n  s  CA(kJ'2Stko2S) 


which  yields  (3.32). 


By  a  similar  procedure  it  is  possible  to  show  also  that  H(R)  >  p 

Remark  3.5.  In  Lemma  3.7,  we  have  assumed  that  s  is  an  integer.  We  can 
generalize  the  results  by  interpolation  of  weighted  spaces,  (see,  e.g.  [7]) 
for  s  nonintegers. 

3.4.  Numerical  examples. 

Consider  now  the  problem  (3.1),  (3.2)  with  £2=1,  f  =  0  and 
2 

i )  gj(x)  =  1  -  x 

ii)  g2(x)  =  1 

iii )  g3(x)  =  cos 

Denoting  u^(t,x)  the  solution  of  our  problem  with  the  initial  function  g^ 
we  can  extend  it  antisymmetrically  with  respect  to  x  =  ±1  into  a  periodic 
function  defined  on  IxlR.  We  will  denote  the  extended  function  once  more  by 

u. .  Let  now  v  =  — .  Then  v  satisfies  the  equation 
3xJ 


av  _  a_v 

at  .2 
ax 


on  Rxl 


v(-l,x) 


-I- 


1)  tf(x-(2k+l) 


where  Q  is  the  Dirac  function.  Hence, 


(x-(2k+l ) ) 


wt.x)  -  C  v  —  (- 


1  ^ 

1 )  e 


and  we  get 


anv 


(t,  x) 


3x 


(l-x2)2+ndx  -  t('n/2)  +  n/2) 


and  hence 


(3.33) 


|u,-pnu,l|2,  <Cp-(3/2+„)2t(-n/2)*(l/2) 

1  0  1  fi'(Q) 


By  interpolation  theory  (3.33)  holds  not  only  for  n  integers  but  also  for  n 
real .  Hence  we  can  choose  n  =  3  -  e  and  get 


(3.34) 

Similarly  we  get 


|ui  “  poui"x  -  c(e)p 


-(4.5-e) 


I  up  -  pnuo  H 0 1  -P 
2  °  2  ft1 (Q) 


.  -(n-l/2)2,(-n/2)-(l/2) 

2  p  t 


and 

(3.35) 


|U2  '  P0U2lfX  -  C(G)P 


(0. 5-e) 


where  in  (3.34)  and  (3.35)  c  >  0  is  arbitrary.  In  case  iii)  we  obviously 
have 

2 


-  *-t 

u3(t,  x)  =  cos  ^e  2 


and  hence 


|u  -  P  u  ||  <  C ( k ) p 

3  °  3 


-k 


for  arbitrary  k  >  0.  Using  the  results  from  [6]  we  get 


W5W 


(3.36) 


|U3  ^^X  “  Cr 


with  r  <  1  and  C  which  can  be  estimated. 

Applying  ow  theorem  3.4  together  with  the  estimates  (3.34),  (3.35),  and 


(3.36)  we  get 


(3.37a) 


|ul ,  s  ~  U1 11 X  5  Cp 


-(4.5-e) 


(3.37b) 


(3.37c) 


1 U2 ,  s  ""  U2 1 X  -  CP 


lu3,s-U3"x"C^ 


-(O.S-c) 


with  c  >  0  arbitrary  and  r  <  1. 

Figure  3.1  shows  the  computed  error  of  ||u.  —  u.  ||  in  the  log  log 

1 ,  S  1  A 

scale.  Figure  3. 1  also  shows  the  slopes  indicated  by  our  estiamtes.  We  can 

see  a  very  good  agreement  for  Uj  and  u^.  In  the  case  u 2  the  numerical 

-1  -1/2 

experiment  seems  to  indicate  a  better  result  (rate  p  and  not  p  )  in 
the  computed  range. 


SPACE  DEGREE  P 


F  i  glare  3.1. 


The  relative  error  for  the  semidiscrete  method  (discretization  in 
x)  with  the  initial  functions: 

2 

(a)  g^(x)  =  1-x  .  slope  based  on  (3.37a); 

(b)  g2(x)  =  sl°Pe  based  on  (3.37b); 

(c)  g3(x)  =  COS  g-. 


3.5.  The  semidiscrete  problem:  Discretization  in  t. 

In  this  section  we  will  consider  the  case  when  we  discretize  the  variable 
t,  while  keeping  x  continuous.  This  is  essentially  the  p-version  of  the 
Rothe  method  (see  [8]).  Let 


T  _j(I)  =  (u  e  L  (I)  |  u  is  a  polynomial  of  degree  q-1) 


.v 


o )  2 

T^(I)  =  {v  €  L  (I)  |  v  is  a  polynomial  of  degree  q  and  v(l)  =  0} 


Define  now 


O  1 


S  =  T  .xH  (n) 
q  q-1 


V  =  f’xftta). 

q  q 


Obviously  c  X  and  c  Y,  where  X  and  Y  are  the  spaces  introduced 


in  section  3.1.  For  any  u  e  S^,  resp.  v  €  V  ,  we  have 


u(t,x)  =  ^«1  (t)u.  (x) 


i=l 


where  u^  are  the  eigenfunctions  as  in  section  3.1,  with 


a.  ( t )  e  T  . 
l  q-1 


and 


=  y  (t)ui (x) 


i  =  l 


with 


IS.  (t)  €  f  ’ . 

i  q 


The  norms  || u|[ x  and  ||v||  are  now  given  by  (3.8)  and  (3.11).  Defining  the 


bilinear  forms  B(u, v)  on  S^xV^  by  (3.13)  we  see  that  (3.14)  holds.  Using 


now  Theorem  2.3  we  get 


inf 

sup  |B(u,v)| 

>  Cq 

veV 

ueS 

q 

q 

-T 

ii 

liu||x<l 

inf 

sup  I B(u, v) | 

>  Cq 

U€S 

veV 

q 

q 

11 

X 

~o 

«— < 

VI 

> 

"> 

•1/2 


■1/2 


41 


,v 

•sa 


si 


SI 


& 


I 


a 


'4 


a 

sa 


% 


:: 


a 


I 


$ 


:3 

& 


Hence  we  have 


Theorem  3.8.  Let  uQ  be  the  solution  of  the  problem  f  and  us  its  semi- 
discrete  solution  (discretization  in  t).  Then 


(3.38) 


lu  -u  ||  <  Cq1/2  inf  ||u-  w|| 

weS 

q 


Theorem  3.9.  Assume  that  f  =  0  in  the  problem  f  and  that  the  conjecture 
A  holds.  Then 


(3.39) 


I us  "  uq II x  ~  c<p(q)  inf  |u-w| 
weS 

q 


where  ip  =  1  or  <p(q)  =  log(q). 

Let  us  remark  that  there  are  solutions  of  the  problem  T  with  f  being 

1/2 

rough  so  that  the  factor  q  in  (3.38)  is  present.  On  the  contrast  if  f 
is  smooth  then  (3.39)  is  applicable.  Let  us  now  consider  the  examples  dis¬ 
cussed  in  section  3.4.  Obviously  we  can  write 


where 


and 


u(t,x)  = 

00 

Zv 

,2k-l . 2  2,.  ,  . 

( — 2~)  71  (t  +  1) 

cos (2k- 

k=l 

lak' 

<  ck-3 

for 

the.  function  u^ 

lak' 

<  ck  1 

for 

the  function  u^ 

a^  =  0  k  =  2,3,...  for  the  function  u^. 
Using  now  (2.28)  we  get 


42 


•  ii  n  2  \  ,2  2 ,2n~ 2  —2n 

inf  || u  -  w||  <  >  k  a  (kit)  q 

W€S  H 

q  k=l 


and  hence  we  can  use  for  u^,  n  =  3-e  and  get 


(3.40a) 


analogously 


(3.40b) 


'“l.s'Vx  S  C<I 


l^.s-Vx  £  C<< 


-(2 -c) 


...  -2.5+e  -0.5+e  .„  .  ,  ,  ,  , 

or  possibly  q  and  q  if  conjecture  A  holds  and 


(3. 40c) 


for  k  arbitrarily  large. 


|U3,s-u3"x  £  Cc! 


Figure  3.2  presents  the  computed  error  of  ||u.  —  u .  ||  in  the  log  log 

1  t  S  1  A 

scale.  Figure  3.2  also  shows  the  slope  indicated  by  our  estimates  taking  into 
account  our  conjecture  A.  We  can  see  a  very  good  agreement  for  all  cases. 


'■V.V.V.V.V.V.V.V.V.V.V.V 


TIME  DEGREE  q 


Figure  3.2.  The  relative  error  for  the  semidiscrete  method  (discretization  in 
t)  with  the  initial  functions: 

2 

(a)  g^x)  =  1-x  ,  slope  based  on  (3.40a)  with  conjecture  A; 

(b)  ~  1  *  slope  based  on  (3.40b)  with  conjecture  A; 

(c)  g3(x)  =  cos 


3.5.  The  complete  discretization. 

In  this  section  we  will  consider  the  general  case  when  the  discretization 
in  x  and  t  has  been  made  simultaneously.  Let  us  define 

S  =  S(R,q)  =  T  xR 


V  =  V(R,q)  =  f’xR. 


Now  we  obviusly  have  for  u  e  S,  v  e  V 


(3. 41a) 


I  B(  u,  v)  |  <  II u|| x(|  v||  y 


said  combining  the  results  of  previous  sections  we  get 

-1/2 


(3.41b) 

inf 

sup 

ueS 

veV 

H 

VI 

X 

"3 

H 

ii 

> 

> 

(3.41c) 

inf 

sup 

veV 

ueS 

II  v||  y=l 

H 

VI 

X 

~3 

where  C 

is  independent  of 

R  and 

-1/2 


Theorem  3. 10.  Let  u^  be  the  solution  of  the  problem  9  and  u^  e  S  be  the 


finite  element  solution 


B(u  , v )  =  F(v)  V  v  e  V. 


Then 


lu  -u0Hx  <  Cq1/2(K(R})_1  inf  ||u  -w|| 

weS 


In  the  case  that  f  =  0  we  get 


I ~  Uq II x  -  C<p(q)  (K(R) )  1  inf  ||u  -w||s 

weS 


where  <p(q)  <  1  or  <p( q)  <  log(q)  when  the  conjecture  A  holds. 


Let  us  now  proceed  for  another  estimate.  Denote  solution  when  the 


variable  t  is  discretized,  while  x  is  continuous.  Let  ||p  -u  ||  <  e  . 

1  U  X  1 


Proceeding  as  before  we  get  now 


|US  "  “o"x  s  c 


[CW  ,"'T'Po',T»gl(n)dt] 


and 


45 


V; 


v1 


a 


m 


Cv 

w 


W 


& 


I 


a 


i 


•  v‘ 


I 


,4 


jn'jnWTKWC  V’JC'3tTTKV«?WJ’V CVJV. -WW.YJWTO 


^T  ~  P0PT  *  ^  U0^  uo^  +  u0  "  P0u0' 


Using  now  Theorem  3.6,  we  get 
.1 


HHj.-PqPUj.I^  «  iC  ||u0-P0u0l|2dt  +  e^(l+7,(R)+X(R))?(R)J 


and  hence 


lus-u0ilx  *  C(c  (l+T»(R)+A(R))€(R)+e2) 


where  e^  is  the  error  of  the  semidiscrete  method  (the  discretization  of  t), 
which  was  discussed  in  section  3.5  and  e^  is  the  error  of  the  semidiscrete 
method  (the  discretization  in  x)  which  was  discussed  in  section  3.2. 

Let  us  consider  now  the  numerical  example  we  discussed  earlier.  In  this 
case  we  have 

(1+t}(R)+A(R))C(R)  *  Cp 

1/2 

but  using  Lemma  3.7  we  can  replace  this  term  by  Cp  .  hence  we  get 

,,  ,  ,,  .,  .  „  1/2-  -5+e  -2. 5+e. 

(3.42a)  u,  -  u,  JL,  <  Cp  (p  +q  ) 


(3.42a) 

n  I,  ^  -  1/2-  -5+e 

l|ul'ul,sl!x  *  Cp  (p  +q 

(3.42b) 

|u2"u2,sBX  -  Cp  (P  +q 

(3.42c) 

1  u3  _  U3 , s 11 X  -  CP1/2(P_k+c! 

where  e  >  0, 

k  <  0  anbitrany 

The  rates  in  q  are  increased  by  1/2  since  the  conjecture  A  has  been 

used.  The  estimates  (3.42a,b,c)  indicate  the  optimal  combination  in  q  and 

p.  Taking  into  account  the  fact  that  the  computational  work  (i.e.,  number  of 

3 

arithmetic  operations)  is  proportional  to  qp  ,  we  can  see  that  the  optimal 
2 

choices  are  q  ~  p  in  the  first  and  second  cases,  and  q  ~  p  in  the  last 


I 

I 


case.  Using  the  estimate  (3.42a,b,c)  we  get  the  optimal  relationship  between 


the  error  and  the  work: 

(3.43a)  E  =  CW~°‘9  for  g(x)  =  1  -  x2 

(3.43b)  E  =  CW-0’1  for  g(x)  =  1 

-If  iry 

(3.43c)  E  =  CW  ,  k  >  0  arbitrary  for  g(x)  =  cos 

Figures  3.3  -  3.5  present  the  performance  of  the  complete  discretization  via 

p-version  of  the  finite  element  method.  Figures  show  the  slopes  indicated  by 

the  functions  given  by  (3.43a)  and  (3.43b).  We  can  see  a  very  good  agreement 

for  u^  (Figure  3.4,  numerical  rate  of  convergence  is  slightly  better  than 

the  theoretical  one)  and  for  u^  (Figure  3.5).  In  the  case  (a),  i.e.,  for 
2 

g^(x)  =  1-x  (Figure  3.3)  the  numerical  rate  of  convergence  is  very  close 

to  the  theoretical  one,  but  the  best  result  is  obtained  for  the  combination 

18  2 
p  =  [q  ]  instead  of  p  =  q  . 
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Figure  3.3.  The  performance  of  the  complete  discretization  via  p-version  of 

the  finite  element  method  for  various  combinations  between  the 

space  degree  p  and  the  time  degree  q.  The  initial  function  is 
2 

g^(x)  =  1  -  x  .  The  slope  is  based  on  (3.43a). 
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Figure  3.4.  The  performance  of  the  complete  discretization  via  p-version  of 
the  finite  element  method  for  various  combinations  between  the 
p  and  q.  The  initial  function  is  g2(x)  =  1  ’  The  sloPe  is 
based  on  (3. 43b) . 
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Figure  3.5.  The  performance  of  the  complete  discretization  via  p-version  of 
the  finite  element  method  for  various  combinations  between  p 

TTX 

and  q.  The  initial  function  is  g^(x)  =  cos—. 


3.7.  Additional  comments. 

In  section  2.6  we  mentioned  that  the  theory  can  be  extended  where  compo¬ 
site  intervals  are  used.  This,  of  course,  can  be  applied  here  too.  For 
example  the  estimate  (3.42a)  indicates  that  the  refinement  in  the  variable  t 
would  be  advantageous. 

We  have  assumed  that  the  discretization  in  x  is  the  same  for  all 


tel.  If  I  is  partitioned  into  I .  then  also  discretization  in  x  can  be 

J 


different  for  various  I Then  we  get  the  h-p  version  for  parabolic  equa¬ 


tion.  These  and  analogous  questions  will  be  discussed  in  forthcoming  papers. 


v  www  j.  w  w  ^ 


til 

t2] 

[3] 

[4] 

lb] 

[6] 

[7] 

[8] 


uc  f  w  vj  toj 


References 


Mackerle,  J. ,  MAKABASE,  An  information  system  on  structure  mechanics 
software  and  applications,  Adv.  Eng.  Software  8  (1986),  81-87. 

Mackerle,  J. ,  MAKABASE,  An  On-line  Information  Retrieval  System  of  Struc¬ 
tural  Mechanics,  Computers  and  Structures  24  (1986),  977-983. 

Babuska,  I.,  The  p  and  h-p  versions  of  the  finite  element  method,  The 
State  of  the  Art  Proc.  of  the  Workshop  on  Finite  Element  Method,  NASA 
Langley  and  R.  Voigt,  Springer  1988,  199-239. 

Babuska,  I.,  Aziz,  A. K. ,  Survey  Lectures  on  the  Mathematical  Foundations 
of  the  Finite  Element  Method  in  The  Mathematical  Foundations  of  the 
Finite  Element  Method  with  Applications  to  Partial  Differential 
Equat ions,  ed.  A. K.  Aziz,  Academic  Press,  New  York,  London,  1972,  3-363. 

Arnold,  D. N. ,  Babuska,  I.,  Osborn,  J. ,  Finite  Element  Methods,  Principles 
for  Their  Selection,  Comp.  Math.  Appl.  Mech.  _  45  (1984),  57-96. 


Gui,  W. ,  Babuska,  I.,  The  h,  p  and  h-p  version  of  the  finite  element 
method  in  1  dimension,  Part  I,  Part  II,  Numer.  Math.  49  (1986), 
577-612,  613-657. 

Triebel,  H. ,  Interpolation  Theory.  Funct ions  Spaces.  Differential 
Operators.  North-Hoi  land,  Amsterdam,  New  York,  Oxford,  1978. 

Rektorys,  K. ,  The  Method  of  Discretization  in  Time.  D.  Reidel  Publishing 
Co. ,  Dordecht,  Boston,  London,  1982. 


> 


50 


A 


•o 


» 


39 


3 


S3 


N 

..V 


A 


The  Laboratory  for  Numerical  analysis  is  an  integral  part  of  the 
Institute  for  Physical  Science  and  Technology  of  the  University  of  Maryland, 
under  the  general  administration  of  the  Director,  Institute  for  Physical 
Science  and  Technology.  It  has  the  following  goals: 

o  To  conduct  research  in  the  mathematical  theory  and  computational 

implementation  of  numerical  analysis  and  related  topics,  with  emphasis 
on  the  numerical  treatment  of  linear  and  nonlinear  differential  equa¬ 
tions  and  problems  in  linear  and  nonlinear  algebra. 

o  To  help  bridge  gaps  between  computational  directions  in  engineering, 
physics,  etc.,  and  those  in  the  mathematical  community. 

o  To  provide  a  limited  consulting  service  in  all  areas  of  numerical 
mathematics  to  the  University  as  a  whole,  and  also  to  government 
agencies  and  industries  in  the  State  of  Maryland  and  the  Washington 
Metropolitan  area. 

o  To  assist  with  the  education  of  numerical  analysts,  especially  at  the 
postdoctoral  level,  in  conjunction  with  the  Interdisoipl inary  Applied 
Mathematics  Program  and  the  programs  of  the  Mathematics  and  Computer 
Science  Departments.  This  includes  active  collaboration  with  govern¬ 
ment  agencies  such  as  the  National  Bureau  of  Standards. 

o  To  be  an  international  center  of  study  and  research  for  foreign 

students  in  numerical  mathematics  who  are  supported  by  foreign  govern¬ 
ments  or  exchange  agencies  (Fulbright,  etc.) 

Further  information  may  be  obtained  from  Professor  I.  Babuska,  Chairman, 
Laboratory  for  Numerical  Analysis,  Institute  for  Physical  Science  and 
Technology,  University  of  Maryland,  College  Park,  Maryland  207*42. 
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